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Abstract 

We present a new measurement of the top quark mass in tt events 
in which both W bosons from top quarks decay into leptons (eu,fj,v). 
We use events collected by the CDF experiment from pp collisions at 
v /s=1.8 TeV at the Tevatron collider. We measure a top quark mass of 
167.4 ± 10.3(stat) ± 4.8(syst) GeV/c 2 from a sample of eight events. We 
combine this result with previous CDF measurements in other decay chan- 
nels to obtain a mass value of 176.0 ± 6.5 GeV/c 2 . 
PACS numbers: 14.65.Ha, 13.85.Qk, 13.85.Ni 
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A precise measurement of the top quark mass is an important ingredient in testing the 
consistency of the standard model with experimental data. In addition, precise W and top 
mass measurements can provide information on the mass of the Higgs boson, which is a 
remnant of the mechanism that gives rise to spontaneous electroweak symmetry breaking. 
Using 109 pb _1 of data accumulated by the CDF experiment at the Fermilab Tevatron from 
1992 through 1995, we report an improved measurement of the top quark mass using dilepton 
events originating predominantly from tt — > W + bW~b — > (£ + vb) (£~ i?b) , where £ = e or \i. 
This measurement supersedes our previously reported result in the dilepton channel 0. The 
previous result was obtained by comparing data with Monte Carlo simulation of tt events 
for two kinematic variables, the 6-jet energies and the invariant masses of the lepton and 
6-jet systems. Here we make use of all available information in the event and obtain a more 
precise measurement. We combine the result from the dilepton channel with those from the 
lepton plus jets channel (tt — > W + bW~b — > (£ + ub)(qq'b)) || and the all-hadronic channel 
(tt — > (qq'b)(qq'b)) |J to obtain an overall top quark mass value from the CDF data. 

The CDF detector consists of a magnetic spectrometer surrounded by calorimeters and 
muon chambers ||. A four-layer silicon vertex detector (SVX) |5[], located immediately 
outside the beampipe, is surrounded by the central tracking chamber (CTC) which is inside 
a 1.4 T superconducting solenoid. This tracking system is used to measure the momenta of 
charged particles. Electromagnetic and hadronic calorimeters, located outside the CTC, are 
segmented in projective towers and cover the pseudorapidity region < 4.2 0. They are 
used to identify electron and photon candidates and jets, and are used to measure the missing 
transverse energy ($t) which can indicate the presence of energetic neutrinos. Outside the 
calorimeters, drift chambers in the region < 1.0 provide for muon identification. A 
three-level trigger selects events that contain a high P? electron or muon for this analysis. 

We apply the same event selection criteria and use identical background calculations as 
those employed in the previous mass analyses of the dilepton channel |1[]. We require two 
high transverse momentum (Pp > 20 GeV/c) oppositely charged leptons (e or /i) in the 
central detector region (|?7| < 1), with at least one of them well isolated from nearby tracks 
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and calorimeter activity. To reject Z — > £ + £~X events we require that the dilepton invariant 
mass, M ee or M M(1 , be outside the interval 75-105 GeV/c 2 , and remove events containing an 
isolated photon with Et > 10 GeV if they are consistent with radiative Z decays. We 
require at least two jets in the region \rj\ < 2.0, each with observed Et > 10 GeV in a cone 
radius AR = y/ Ar/' 2 + A</> 2 = 0.4, and require \$t\ > 25 GeV as a signature for missing 
neutrinos. To reject events in which $t is due to lepton or jet energy mismeasurements we 
require \Et\ > 50 GeV if $t is close to a lepton or a jet (A0(E T J or j) < 20°). Finally, 
we require Ht > 170 GeV, where Ht is the scalar sum of the \Pt\ of the two leptons {Et 
for electrons), the Et of the two highest Et jets and \Et\- We obtain a sample of eight 
candidate events. The expected background of 1.3 ± 0.3 events consists of events in which a 
track or a jet is misidentified as a lepton (0.29 events), Drell-Yan production (0.35 events), 
WW production (0.24 events), Z — > rr decays (0.26 events) and Z — > decays in which 
(jl tracks are mismeasured (0.20 events). 

To extract a top mass measurement from these events the jet energies are corrected for 
losses in cracks between detector components, absolute energy scale, contributions from the 
underlying event and multiple interactions, and losses outside the clustering cone. These 
corrections are determined from a combination of Monte Carlo simulation and data . An 
additional energy correction is applied to jets containing a muon to take into account the 
low calorimeter response to muons (we add the Et of the muon to the energy) and missing 
neutrinos from semileptonic decays of heavy quarks in the jets. We obtain a corrected Bt 
by correcting each of its components, i.e., the jets, the leptons and the energy which has 
not been clustered into jets. 

The procedure used here has two steps: we reconstruct each event to obtain a top mass 
estimate for the event, then we apply a likelihood method to obtain an overall mass value 
in the presence of background. Each candidate event is reconstructed according to the ti 
decay hypothesis in the dilepton channel: 

t -> W + b -> i\v x b 
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t -> W~b -> £~z/ 2 6 

The two highest jets in the event are assumed to be the 6-jets from top decays. We assume 
the 6-jet mass to be 5 GeV/c 2 . After applying the invariant mass constraints m(i\V\) = 
^(^2^2) = triw and m{£\V\b) = m^^b), the system remains underconstrained due to 
the two unmeasured neutrinos. Therefore, for any assumed top mass value mt, we use 
a weighting technique to determine a function, f(m t ), from which we extract a top mass 
value 0. We proceed as follows. We assume a top quark mass (m t ) and the two neutrino 77 
values (771, rj2, see below) and solve for the neutrino momenta, up to a four- fold ambiguity 
(two P z {y) choices for each v) for each of the two jet charged- lepton pairings. We then assign 
a weight to each solution by comparing Et P , the sum of the neutrino transverse momenta 
for that solution, to Et"" , the measured missing transverse energy after proper correction: 

g[rnt,T}uV2) = exp\ — I 

where a is the Et resolution for that event (see below). 

For each choice of m t , 771, and 772 we take into account the detector resolution for jets 
and leptons by sampling (i.e. fluctuating) the measured quantities many times accord- 
ing to their resolutions. The Gaussian resolutions for electrons and muons are oeIE = 
^/(0.135) 2 /£ T + (0.02) 2 and a(l/P T ) = 0.11%, where E T , in GeV, is the electron energy 
measured in the calorimeter and Pt is the beam-constrained muon momentum measured in 
the CTC. For jets we use an E T dependent Gaussian resolution appropriate for 6-partons 
derived from the HERWIG ti Monte Carlo, in conjunction with detector simulation, as- 
suming a top mass of 170 GeV/c 2 0. The Et is recomputed for each sampling using the 
new jet and lepton energies. This procedure takes into account all the uncertainties in the 
Et measurement except for the resolution of the unclustered energy measurement, which 
is estimated to be 4 GeV for both transverse components for low luminosity events. Hence 
we use a= 4-y/rT GeV in the above expression for g(m t , r/ 1? r/ 2 ), where n is the number of 
interactions in the event. This value has been obtained from minimum bias events and for 
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each event is properly scaled to take into account the effect of multiple interactions at high 
luminosity. 

For each assumed top mass value we use several (100) values of the two neutrino rji 
and 7/2, chosen from distributions obtained from the HERWIG Monte Carlo predictions ||. 
They are consistent with independent Gaussian distributions with a = 1.0 unit of r/. The 
weight is summed for all samplings as well as over all rji, 772 values and all the eight possible 
combinations; thus for each event at each top mass, m t , we evaluate an overall weight: 

f( m t) = H gimtiVuVz) 

where refers to all the jets in the event. We then compute the weight as a function of 
the top mass in 2.5 GeV/c 2 steps in the range 90-290 GeV/c 2 . 

The f{m t ) distribution for each of the eight candidate events, normalized to unity, is 
shown in Figure [I]. For each event, i, we use this distribution to determine a top mass 
estimate, rrii, by averaging the values of rrit corresponding to values of f{mt) closest to and 
greater than f{m t ) max /2 on either side of the maximum. The rrii distribution for the eight 
events is shown in Figure |2|, together with the Monte Carlo expectation for background 
alone, and top plus background normalized to the data. 

Given the rrii distribution for a sample of N events, we use a maximum likelihood method 
similar to that employed in the lepton plus jets mass analysis || to extract a top mass value. 
We write the likelihood as: 

r n< - \ v(at 1 \ tt n s T{rrii, m t ) + n b B(mi) 

L = G(n b ; n b , a b ) ■ P(N; n s + n b ) ■ M 

\ n s + n b 

where T(rrii, m t ) and B(mj) are the probability density functions (templates) for recon- 
structing a mass rrii from tt events with a true top mass m t and from background events, 
respectively; n s and n b are the numbers of tt and background events; G(n b ; n b , a b ) is a Gaus- 
sian in n b with mean n b (1.3 events) and width o b (0.3 events) and P(N; n s + n b ) is a Poisson 
distribution in N with mean n s + n b . The top templates, T(rrii,m t ), are obtained by per- 
forming the dilepton mass reconstruction algorithm on HERWIG tt Monte Carlo samples, 
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and parameterizing the resulting reconstructed mass distributions as smooth functions of 
both rrii and m t 0. The background templates, B(rrii), are derived from data and Monte 
Carlo ]!]. The likelihood fit gives the top mass measurement and a statistical error taken 
from the change in the -log(£) of 0.5 units from its minimum value. 

To check the likelihood fitting technique, we use Monte Carlo events to perform many 
pseudo-experiments for several input top mass values. At each top mass value each experi- 
ment consists of a total of eight top and background events, with the number of background 
events drawn from a binomial distribution with a mean of 1.3±0.3 events. The mass estimate 
for each event is drawn from the appropriate parameterized template. For each experiment 
we obtain a mass, m exp , and a statistical error, a exp . We study the distributions of the pulls, 
{ m exp — m t)/cexp, and find that the medians are consistent with zero and the widths are at 
most 1.1. We take this width into account by multiplying the error returned by the fit by 
a factor 1.1. This factor is included in all statistical errors on the dilepton measurements 
given below. 

We apply the likelihood method to the data shown in Figure |2|. The inset shows the 
negative log-likelihood as a function of the top mass, from which we determine a top quark 
mass value of 167.4 ± 10.3 (stat) GeV/c 2 . Monte Carlo studies performed with pseudo- 
experiments yield an 8% probability for one such experiment to have a statistical error 
< 10.3 GeV/c 2 at a mass of 167 GeV/c 2 . In comparison, the previous result was 161 ± 17 
(stat) GeV/c 2 @. 

The mass reconstruction procedure used here is quite different from that used for the 
lepton plus jets sample @. In that case there is only one missing neutrino, and a kinematic 
fit with two constraints can be performed. We cross-check the present procedure by applying 
it to the five events of the lepton plus four jets mass sample in which two jets are tagged as b- 
jets by the SVX . We assume that the two untagged jets are the products of the hadronic 
W decay and mimic the dilepton decay by treating the highest Et untagged jet as a lepton 
and the second untagged jet as a neutrino. The present likelihood method then gives a top 
quark mass from the five events of 181.5 ± 12.6 GeV/c 2 , which differs by 11 GeV/c 2 from 
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that obtained with the lepton plus jets kinematic fit procedures, 170.1 ± 9.3 GeV/c 2 0. In 
comparison, Monte Carlo studies show that the difference in mass obtained with the two 
methods is expected to be centered at zero with a resolution of 14 GeV/c 2 . 

The systematic errors on the mass measurement are estimated with the same general 
procedure as in the lepton plus jets mass analysis 0: we generate new top mass distributions 
by varying the appropriate quantities in the Monte Carlo simulation and then performing 
likelihood fits to many experiments of eight events each using the standard templates. The 
mass shifts obtained determine the systematic errors, which are summarized in Table |. The 
largest contribution comes from the systematic uncertainties on the jet energy. The second 
largest source, the uncertainty in modeling initial or final state radiation in ti events, is 
estimated using events generated with the PYTHIA Monte Carlo program [[Kj] to isolate 



the effects on the top mass due to initial and final state radiation jets. The error due 
to background shape is estimated by refitting the data with three different backgrounds: 
W + W~, fake leptons, and the sum of the two with the appropriate weights. That due to 
the choice of parton distribution functions (0.6 GeV/c 2 ) is estimated by using an alternative 
set of functions in the HERWIG Monte Carlo. The effect due to the choice of Monte Carlo 
generator (0.6 GeV/c 2 ) is estimated by comparing the mass value obtained with events from 
the HERWIG and PYTHIA Monte Carlo samples. The total systematic error amounts to 
4.8 GeV/c 2 . 

Table | includes the systematic errors for the other two decay channels. The system- 
atic error due to hard gluon radiation uncertainty in the lepton plus jets and the all- 
hadronic || channels are now evaluated in the same way as described above. For the 
all-hadronic case the old value of 8.0 GeV/c 2 was an overly conservative estimate. Also, the 
"fit procedure" systematic in the all-hadronic channel (5.0 GeV/c 2 ) has been removed as 
being overly conservative. The overall systematic error for the lepton plus jets channel has 
increased from 4.9 to 5.3 GeV/c 2 , and that for the all-hadronic channel is reduced from 12.0 
to 5.7 GeV/c 2 . The new top mass measurement for the dilepton channel and those for the 
other two channels with revised systematic errors are shown in Table ||. 
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The results for the three channels are combined with standard methods (TT| to yield an 
overall CDF mass measurement. The three statistical errors are taken as uncorrelated, while 
the systematic errors are assumed to be either entirely correlated or uncorrelated between 
any two channels. The primary systematic error, that due to jet energy uncertainty, is taken 
as entirely correlated among all channels, as is the systematic error due to the Monte Carlo 
model used (gluon radiation and simulation entries in Table |). The combined result is 

m t = 176.0 ± 6.5 GeV/c 2 

including both statistical and systematic errors. In Table we show the combined value 
with separate statistical and systematic errors. They are obtained by defining the combined 
statistical error as the sum in quadrature of the weighted individual statistical errors, and 
the systematic error as the difference in quadrature of the total and statistical errors. The 
relative contributions from the three channels are 67% for lepton plus jets, 18% for dileptons 
and 15% for all-hadronic. 
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TABLES 

TABLE I. Systematic errors (in Gev/c 2 ) on the top mass measurement for the three tt decay 
channels. The Monte Carlo modeling term includes effects from: parton density functions and 
b-tag uncertainty where applicable. For the lepton plus jet channel the Monte Carlo statistics 
term was included in the statistical error. 



Channel 




dilepton 


£+jets 


all-had. 


Jet energy scale 




3.8 


4.4 


5.0 


Initial and final state radiation 


2.7 


2.6 


1.8 


Monte Carlo modeling 




0.6 


0.5 


0.2 


lVTontp (~!arlo p'enprator 




0.6 


0.1 


0.8 


Background shape 




0.3 


1.3 


1.7 


Monte Carlo statistics 




0.7 


n.a. 


0.6 


Total 




4.8 


5.3 


5.7 


TABLE II. 


Summary of top mass measurements with the CDF detector 


Channel 




Top mass (GeV/c 2 ) 




Reference 


Dilepton 




167.4 ± 10.3 ± 4.8 




this paper 


Lepton + jets 




175.9 ± 4.8 ± 5.3 


I 


|], this paper 


All-hadronic 




186.0 ± 10.0 ± 5.7 


1 


3], this paper 


Combined 




176.0 ± 4.0 ± 5.1 




this paper 
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FIG. 1. Weight distribution normalized to unity as a function of m t for the eight dilepton top 
candidate events (a-h). 
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FIG. 2. Reconstructed top mass for the eight dilepton events (solid). Background distribution 
(shaded, 1.3 events) and top Monte Carlo (6.7 events) added to background (dashed). The negative 
log- likelihood distribution as a function of the top mass is shown in the inset. 
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